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Magnetic  nanofluids  (MNF)  constitute  a  special  class  of  nanofluids  that  exhibit  both  magnetic  and  fluid 
properties.  The  interests  in  the  use  of  MNF  as  a  heat  transfer  medium  stem  from  a  possibility  of 
controlling  its  flow  and  heat  transfer  process  via  an  external  magnetic  field.  This  review  presents  recent 
developments  in  this  field  with  the  aim  of  identifying  major  affecting  parameters  and  some  novel 
applications.  This  review  emphasizes  on  thermal  conductivity  enhancement  and  thermomagnetic 
convection  in  devices  using  MNFs  as  heat  transfer  media. 
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1.  Introduction 

Dispersion  of  nano-sized  particles  of  different  materials 
(metals,  metal  oxides,  etc.)  in  a  carrier  fluid  known  as  nanofluids 
has  been  a  subject  of  intensive  investigations  over  decades  due  to 
their  potential  applications  in  heat  transfer  and  electronic  cooling 
[1-6],  Magnetic  nanofluids  (or  ferrofluids),  which  consist  of 
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colloidal  mixtures  of  superparamagnetic  nanoparticles  suspended 
in  a  nonmagnetic  carrier  fluid,  constitute  a  special  class  of 
nanofluids  that  exhibit  both  magnetic  and  fluid  properties 
[7-9],  To  prevent  the  aggregation  (due  to  London-van  der  Waals 
interaction  and  magnetic  interaction  between  the  particles),  the 
suspended  nanoparticles  are  coated  by  a  surfactant  layer  such  as 
oleic  acid  [10-12],  tetramethylammonium  hydroxide  [13],  etc. 
A  wide  range  of  carrier  fluids  are  used,  and  some  magnetic 
nanofluids  are  commercially  available  to  satisfy  different  applica¬ 
tions  [14],  Theoretically,  it  should  be  possible  to  produce  disper¬ 
sion  in  any  liquid  thereby  being  able  to  tailor  the  requirements  of 
viscosity,  surface  tension,  temperature  and  oxidative  stability, 
vapor  pressure,  stability  in  hostile  environments  [15].  However, 
the  choice  of  a  carrier  fluid  for  MNF  suitable  for  heat  transfer 
applications  needs  some  additional  requirements  such  as  high 
conductivity,  high  heat  capacity,  high  thermal  expansion  coeffi¬ 
cient,  etc.  Conventional  heat  transfer  fluid  (such  as  water,  oils, 
ethylene  glycol,  etc.)  could  be  a  superior  option  for  advanced 
applications.  Magnetic  nanoparticles  used  in  magnetic  nanofluids 
are  usually  prepared  in  different  sizes  and  morphologies  from 
metal  materials  (ferromagnetic  materials)  such  as  iron,  cobalt, 
nickel  as  well  as  their  oxides  (ferrimagnetic  materials)  such  as 
magnetite  (Fe304),  spinel-type  ferrites,  etc. 

MNF  has  been  used  as  an  advanced  functional  material  to  position 
the  colloids  at  a  certain  part  of  devices  by  means  of  magnetic  forces. 
For  this  reason,  the  main  aspects  for  research  involve  the  increase  of 
fluid  magnetization,  the  specific  properties  of  free  magnetic  nanofluid 
surfaces  in  the  presence  of  an  external  magnetic  field  and  the 
magnetoviscous  and  magnetorheological  effects  [15].  The  use  of 
MNF  as  transfer  media  becomes  of  importance  in  applications 
[16-18].  In  particular,  a  possibility  to  induce  and  control  the  heat 
transfer  process  and  fluid  flow  by  means  of  an  external  magnetic  field 
opened  a  window  to  a  spectrum  of  promising  applications  including 
magnetically  controlled  thermosyphons  for  technological  purposes, 
enhancement  of  heat  transfer  for  cooling  of  high  power  electric 
transformers,  and  magnetically  controlled  heat  transfer  in  energy 
conversion  systems  [15].  The  heat  transfer  enhancement  using 
MNF  in  the  presence  of  magnetic  field  can  be  classified  as  a 
compound  heat  transfer  technique  with  the  additives  (i.e.,  magnetic 
nanoparticles)  and  an  external  magnetic  field  to  increase  the  heat 
transfer  process  [19].  Compared  to  the  conventional  nanofluids 
(nonmagnetic  nanofluids),  the  use  of  MNF  affected  by  external 
applied  magnetic  field  for  heat  transfer  enhancement  offers  the 
following  advantages  [20]: 

a)  The  possible  absence  of  any  moving  parts  necessary  for 
making  the  fluid  to  flow  in  commonly  energy  conversion  and 
cooling  devices.  The  current  of  MNF  is  generated  by  tempera¬ 
ture  difference  and  non-uniform  magnetic  field,  which  can  be 
formed  by  means  of  a  permanent  magnet  system.  The  config¬ 
uration  of  this  system  determines  the  direction  and  the  type  of 
the  fluid  flow.  As  a  result,  the  thermomagnetic  convection  is 
readily  handled; 

b)  The  thermomagnetic  convection  is  much  more  intensive  than 
the  gravitation  one; 

c)  The  possibility  of  tuning  thermophysical  properties  (thermal 
conductivity  and  viscosity)  of  MNF  using  external  magnetic 
fields  [21,22], 

2.  Preparation  of  MNF 

MNF  is  prepared  via  the  dispersion  of  nano-sized  super- 
paramagnetic  particles  into  a  nonmagnetic  carrier  fluid  such  as 
water,  ethylene  glycol,  hydrocarbon  oil,  etc.  [23],  MNF  used  in 
heat  transfer  applications  is  subjected  to  a  magnetic  field, 
magnetic  field  gradient  and/or  gravitational  field,  which  may 


contribute  to  the  particle  sedimentation  in  the  fluid.  Since  the 
interaction  range  of  magnetic  nanoparticles  in  the  applied  fields  is 
directly  related  to  the  particle  size/size  distribution  of  magnetic 
nanoparticles  [21],  it  is  obvious  that  the  later  plays  a  vital  role  in 
the  particle  sedimentation,  thus  affecting  the  stability  of  MNF. 
The  stability  against  the  particle  sedimentation  may  be  ensured 
when  the  thermal  energy  of  the  particles  becomes  greater  than 
that  of  magnetic  and  gravitational  energies,  respectively,  the 
maximum  particle  size  was  determined  by  Odenbach  [21]  to  be 
d  <  (6kBT//(0M07rH)1/3  for  MNF  used  in  the  presence  of  magnetic 
field  and  d  <  (kBTIApgfm)'13  for  MNF  in  the  presence  of  gravita¬ 
tional  field,  where  kB,  T,  M0,  A p,  g,  d,  H  and  f.i0  denote  the 
Boltzmann  constant,  temperature,  spontaneous  magnetization  of 
the  magnetic  material,  density  difference  between  magnetic 
particle  and  the  carrier  fluid,  gravitational  acceleration,  the  height 
of  the  sample,  the  magnetic  field  and  the  vacuum  permeability, 
respectively.  In  addition,  the  aggregation  of  magnetic  nanoparti¬ 
cles  during  the  synthesis  has  to  be  avoided  at  all  costs.  In 
principle,  the  aggregation  of  particles  increases  their  active 
diameter  and  thus  causes  a  destabilization  of  the  suspension  by 
sedimentation.  The  maximum  particle  i^meter  (d),  in  this  case, 
was  estimated  as  d<  (]44kBT / p0M^J  corresponding  to  the 
maximum  interaction  energy  when  two  interacting  particles 
come  into  contact  [21].  Recent  efforts  have  been  made  to 
synthesize  metal  and  metal  oxide  magnetic  nanoparticles 
with  the  desired  size/size  distribution  [24].  Metallic  nanoparticles 
such  as  Ni,  Fe  and  Co  were  prepared  via  the  techniques  like 
simple  reduction  of  metal-salts,  gas-phase  reduction  of  metal 
complexes,  thermolysis  of  metal-polymer  complexes,  thermal 
decomposition  of  metal-carbonyl  complexes  and  submerged  arc 
nanoparticle  synthesis  system  (SANSS)  [24,25].  Magnetic  nano¬ 
particles  of  metal  oxide  such  as  Fe304,  y-Fe203  and  spinel-type 
ferrites  of  the  formula  MFe204  (with  M  =  Mn,  Co,  Zn,  Ni,  etc.)  are 
mostly  used  in  MNF  due  to  their  chemical  stability.  Metal  oxide 
magnetic  nanoparticles  are  usually  prepared  by  chemical  co¬ 
precipitation,  micro-emulsion  and  recently  phase  transfer 
[14,24,26], 

The  use  of  the  required  superparamagnetic  particles  in  MNF  is 
not  a  definitive  requirement  of  magnetic  nanofluid  stability  [8[.  A 
suspension  with  magnetic  nanoparticles  in  carrier  fluids  will  not 
be  stable  due  to  the  presence  of  London-van  der  Waals  and 
magnetic  forces,  leading  to  the  irreversible  aggregation  of  the 
particles  and  their  subsequent  sedimentation.  Therefore,  the 
preparation  of  stable  MNF  requires  an  introduction  of  repulsive 
forces  between  magnetic  nanoparticles  to  counteract  the  London- 
van  der  Waals  and  the  dipole-dipole  magnetic  interactions.  The 
repulsive  mechanism  between  the  particles  can  be  achieved, 
either  by  coating  the  particles  with  a  polymer  surfactant,  which 
produces  an  entropic  repulsion,  and/or  by  charging  the  surface 
of  the  particles,  producing  a  coulombian  repulsion  [8,14].  It  is 
interesting  to  note  that  the  selection  of  the  mechanism  to  be  used 
should  depend  mainly  on  the  properties  of  the  carrier  fluids  and 
the  particles.  The  dispersion  process  is  usually  performed  in  the 
presence  of  a  polymer  surfactant  by  ultrasonic  equipment  and/or 
a  high  speed  homogenizer. 


3.  Thermal  conductivity  of  MNF 

The  original  idea  of  using  a  suspension  for  heat  transfer 
application  was  a  possibility  of  enhancing  the  thermal  conductiv¬ 
ity  of  common  heat  transfer  fluids  by  the  addition  of  nanoparti¬ 
cles  with  a  higher  thermal  conductivity  [4[.  For  this  reason,  some 
previous  investigations  on  thermal  conductivity  of  nanofluids 
were  dominated  by  nanofluids  prepared  with  metallic  or  metallic 
oxide  nanoparticles  such  as  Ti02,  A1203,  Cu,  CuO,  Ag,  carbon 
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Table  1 

Summary  of  frequently  used  models  for  thermal  conductivity  of  MNF. 


Model  Expression 


Description 


Maxwell-Gannet 
model  [70] 
Modified 
Maxwell- 
Gannet  model 
[69] 

Microconvection 
model  [71] 
Bruggeman 
model  [72] 


Jeffrey  model 
[73] 


Rayleigh  model 

[74] 

Murshed  model 

[75] 


Kp+2Kf+2(p(Kp-Kf ) 
Kp+2Kf-(p(Kp-Kf) 

Kj  _  Kp  +2Kf  +2<pj(Kp—  Kf) 

Kf  ~  Kp+2Kf-(pi(Kp-Kf) 


(1  +ARe2Pr°  ">  [KC^-lZ-v\ 
i  =  l[(39>-l)^+(2-3<p,]+|v/A 

WithA  =  (3<p-l)2(^)2  +  2(2  +  9<p-9<p2)(|) 
i  =  1  +  3<p  +  Cp2  [3k2  +  aji  +  £  (*$)  +  ^  +  ■  •  ] 

With  k  =  and  a  =  ^ 

a  +  2  Kf 

K/  =  Kf  +  3  <p  2K,  +Ktqip  +3939tir‘(KpKl)/4Kf  +  3 KpJK„Kf 


i/  f  V  0-r>"’(KP  -l1  *  l)(Kp  I  2a)K/)y^[Q,y3(aJ- 1)+  1 11 

_  \  I  xyjpA+i)  / 


Qp>'6K/ 


3A2  + 


3A 


9  A  /  Kcp+2Kf 

16  V2Kcp  +  3/fy 


3  A“ 


+  \  yPcp'-p-cpds 


3k6T(l-1.5y3Qp)  |  Gt 
27icpy3r3  +  6nrpds 


Spherical  particles  with  low  volume  fraction  (p 
Spherical  particles  with  volume  fraction^,;  i  =  x,y,z 


A  is  a  constant,  Reand  Prare  Reynolds  and  Prandlt  numbers 
respectively,  and  yis  a  system  dependent  coefficient. 
Spherical  particles,  with  high  particle  volume  fraction  cp 


High  order  terms  represent  pair  interaction  of  randomly  dispersed 
particles 


Suspensions  of  spherical  particles  with  a  regular  particle  distribution 

Murshed  model  includes  the  effects  of  particle  size,  nanolayer,  particle 
movements,  interactions,  and  surface  chemistry  of  nanoparticles 


With  y  =  1  +  3-,  y^=^+tA=l 


T^TZK-f' 


Kcp=K,r  (K'r-KpHy^C+Kp)’  ^cpand  cp~cp  are  density  and  specific  heat  of 
complex  particle,  respectively,  and  they  are  given  by  Pcp  =  ^Pp+ 

3^(2t2+2brpt+b2r^)-3^[(2+2b+b2)+brp(ftrp+2bt+o] 


Cn-cv  — 


^cp-p+(i-^) 


3CJ’=*(2t2+2  bYpt+b'2!-2)-^ 


pb 


pb 


[t2(2+2b'+b'2)b'rp(b'rp+2b't+2t)j 


Cp  pdnd  Cp  | stand  for  the  specific  heat  capacity  at  constant  pressure 
for  nanoparticles  and  carrier  fluid,  respectively.  GTis  the  total  potential 
energy  between  two  interacting  colloidal  nanoparticles. 


wherep  =  3 r2 + 3 rp t  + 12 ; b  =  In b'  =  In ) 


nanotubes  [4,27-30].  These  affected  investigations  on  thermal 
conductivity  of  MNF  since  common  magnetic  materials  used  in 
MNF  have  a  relatively  low  thermal  conductivity,  especially  in  the 
absence  of  external  magnetic  field.  However,  recent  investiga¬ 
tions  have  proven  that  a  solid  material  with  a  high  thermal 
conductivity  is  not  always  an  effective  option  for  a  suspension 
to  enhance  the  thermal  conductivity  of  a  carrier  fluid  [31].  This 
finding  as  well  as  a  possibility  of  controlling  the  thermal  con¬ 
ductivity  of  MNF  via  an  external  magnetic  field  could  be  plausible 
reasons  of  recent  great  interests  in  the  enhanced  thermal  con¬ 
ductivity  of  MNF.  The  thermal  conductivity  enhancement  is 
characterized  by  a  thermal  conductivity  ratio,  which  is  defined 
as  a  quotient  of  the  thermal  conductivity  of  the  (magnetic) 
nanofluid  (K)  to  that  of  the  carrier  fluid  (Kf). 

3.3.  Experimental  investigations  on  thermal  conductivity 

Experimental  investigations  on  thermal  conductivity  of  MNF 
are  mostly  conducted  at  room  temperature  with  the  common 
methods  for  conventional  (nonmagnetic)  nanofluids  [32].  The 
mostly  used  techniques  for  MNF  include  the  transient  hot  wire 
[13,33-35]  and  thermal  constants  analyzer  techniques  [22,36,37]. 
Paul  et  al.  [32]  gave  the  detailed  information  regarding  these  two 
techniques  as  well  as  their  advantages  and  disadvantages. 


3.3.3.  Thermal  conductivity  enhancement  in  the  absence  of  magnetic 
fields 

Experimental  investigations  on  the  thermal  conductivity  of 
MNF  in  the  absence  of  magnetic  fields  show  that  the  thermal 
conductivity  enhancement  is  mainly  affected  by  different  para¬ 
meters,  i.e.,  volume  fraction  of  magnetic  nanoparticles,  particle 
size/particle  size  distribution,  chemical  composition  of  magnetic 
nanoparticles,  temperature,  particle  coating  layer,  etc. 

Some  work  on  the  effect  of  particle  volume  fraction  showed 
that  the  thermal  conductivity  of  MNF  increases  with  the  particle 
volume  fraction  [13,22,33,34,37-40].  Abareshi  et  al.  [13]  mea¬ 
sured  the  thermal  conductivity  of  a  water  based  magnetite 
nanofluid  as  a  function  of  particle  volume  fraction  at  different 
temperatures.  The  thermal  conductivity  increased  with  the 
increase  of  the  particle  volume  fraction  and  temperature.  The 
highest  thermal  conductivity  ratio  observed  was  11.5%  at  a 
particle  volume  fraction  of  3%  at  40  °C.  Their  experimental  data 
showed  a  fairly  good  agreement  with  the  simulated  results  by  the 
Murshed  model  (see  Table  1).  Li  et  al.  investigated  the  effects  of 
particle  volume  fraction,  surfactants  and  magnetic  field  on  the 
transport  properties  of  a  water  based  Fe  magnetic  nanofluid  [33]. 
They  observed  a  thermal  conductivity  enhancement  of  MNF  with 
the  increase  of  particle  volume  fraction  both  with  and  without 
applied  magnetic  fields.  Their  results  also  showed  that  the 
viscosity  of  MNF  increased  with  the  particle  volume  fraction 
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of  the  suspended  magnetic  nanoparticles  and  the  surfactants.  Yu 
et  al.  [38]  investigated  the  effects  of  particle  volume  fraction  on 
the  thermal  conductivity  of  a  kerosene  based  Fe304  magnetic 
nanofluid  prepared  via  a  phase-transfer  method.  The  thermal 
conductivity  ratios  obtained  increased  linearly  with  the  increase 
of  volume  fraction  and  temperature  and  the  value  was  up  to 
34.0%  at  1  vol%.  To  further  investigate  the  influence  of  tempera¬ 
ture  on  thermal  conductivity  enhancement,  the  measurements 
were  performed  in  the  temperature  range  from  10  to  60  °C.  The 
results  revealed  that  the  absolute  thermal  conductivity  increased 
with  increasing  temperature,  while  the  thermal  conductivity  ratio 
was  almost  constant  and  the  thermal  conductivities  of  the  MNF 
tracked  those  of  the  carrier  fluid.  Philip  et  al.  [22,37]  conducted  a 
work  on  thermal  conductivity  of  MNF  prepared  via  the  dispersion 
of  magnetite  nanoparticles  in  water,  ethylene  glycol  and  kerosene 
as  carrier  fluids,  respectively.  Their  results  showed  that  the 
thermal  conductivity  ratio  increased  with  the  increase  of  particle 
volume  fraction,  but  there  was  no  thermal  conductivity  enhance¬ 
ment  for  MNF  at  a  volume  fraction  of  <  1.71  vol%.  The  slope  of  the 
linear  region  was  0.035  and  the  highest  thermal  conductivity 
ratio  observed  was  23%  at  7.8  vol%.  Hong  et  al.  [31]  investigated 
the  thermal  conductivity  enhancement  of  an  ethylene  glycol  Fe 
based  magnetic  nanofluid.  They  found  that  the  thermal  conduc¬ 
tivity  ratio  increased  nonlinearly  with  the  increase  of  volume 
fraction.  Hong  et  al.  [39]  investigated  the  thermal  conductivity  of 
nanofluids  with  different  volume  fractions  of  Fe  nanoparticles  in 
ethylene  glycol.  Their  results  confirmed  the  intensification  of 
thermal  conductivity  with  the  particle  volume  fraction.  In  the 
comparison  of  the  copper  and  iron  nanoparticles  dispersed  in 
ethylene  glycol,  the  thermal  conductivity  enhancement  in  iron- 
based  nanofluids  was  higher  than  that  in  copper-based  nanofluid. 
The  observed  thermal  conductivity  enhancement  was  attributed 
to  the  fact  that  the  nanoparticles  in  the  fluids  with  high  volume 
fractions  formed  clusters  at  a  higher  rate.  Zhu  et  al.  [41]  prepared 
a  distilled  water  based  Fe304  magnetic  nanofluids  and  investi¬ 
gated  the  effect  of  volume  fraction  on  the  thermal  conductivity 
enhancement  of  the  MNF.  Gutierrez  [34]  analyzed  the  thermal 
conductivity  of  water  based  spinel-type  ferrite  magnetic  nano¬ 
fluids  prepared  with  magnetic  nanoparticles  of  MnxZn(1_x)Fe204 
(with  x=0,  0.3,  0.7,  1).  The  analysis  of  the  volume  fraction 
dependence  on  the  thermal  conductivity  showed  that  the  thermal 
conductivity  ratio  increased  linearly  with  the  increase  of  volume 
fraction.  In  addition,  the  investigations  on  the  particle  composi¬ 
tion  effects  revealed  that  the  thermal  conductivity  ratio  does  not 
depend  on  the  chemical  composition  of  magnetic  nanoparticles. 
The  enhancement  of  thermal  conductivity  was  attributed  to  the 
effects  of  the  particle  size  and  the  volume  fraction  of  magnetic 
nanoparticles.  The  analysis  of  the  carrier  fluid  effects  showed  that 
the  thermal  conductivity  ratio  is  higher  for  a  carrier  fluid  with  a 
low  thermal  conductivity  such  as  common  hydrocarbon 
(K=0.12-0.14  W/I<  m).  However,  the  absolute  thermal  conduc¬ 
tivity  of  a  MNF  is  higher  for  a  carrier  fluid  with  a  high  thermal 
conductivity.  For  instance,  the  absolute  thermal  conductivity  for 
water  based  magnetic  nanofluid  (K=0.801  W/Km)  reported  was 
three  times  greater  than  that  of  oil  based  one  (I<=0.255  W/K  m) 
for  the  highest  volumetric  concentration  of  both  fluids.  Holotescu 
et  al.  [42]  reported  on  the  utilization  of  the  semi-empirical 
equation  for  the  effective  thermal  conductivity  of  the  magneto- 
granulometric  analysis  to  estimate  the  thermal  conductivity  of 
transformer  oil  based  magnetite  magnetic  nanofluids.  In  their 
analysis,  the  semi-empirical  equation  for  the  effective  thermal 
conductivity  of  the  Holotescu-Stoian  model  [43]  was  applied, 
using  the  results  obtained  for  the  size  distributions  from  the 
magnetogranulometry  analysis,  followed  by  a  comparison  with 
the  measured  values  of  the  effective  thermal  conductivity 
obtained  by  a  hot  ball  method.  It  was  revealed  that  the  use  of 


the  magnetogranulometric  approximation  of  the  magnetite  nano¬ 
particle  size  distribution  could  evaluate  the  thermal  conductivity 
more  efficiently,  compared  to  the  use  of  the  classical  Maxwell 
equation.  They  also  ascertained  the  importance  of  the  data 
regarding  the  analysis  of  the  lognormal  distribution  by  either 
number  or  volume,  as  well  as  the  fact  that  the  results  were 
strongly  dependent  on  the  experimental  method  and  the 
associated  model. 


3.1.2.  Thermal  conductivity  enhancement  in  the  presence  of 
magnetic  field 

Recent  work  on  the  control  of  the  thermal  conductivity  of  MNF 
has  received  a  particular  interest  because  the  “anomalous" 
thermal  conductivity  observed  under  the  influence  of  an  external 
magnetic  field.  In  the  presence  of  magnetic  fields,  besides  the 
aforementioned  parameters  affecting  the  thermal  conductivity  of 
MNF  without  magnetic  fields,  the  thermal  conductivity  of  MNF 
can  be  affected  by  the  orientation  and  the  intensity  of  the  applied 
magnetic  field  [15,33].  Measurements  of  thermal  conductivity  in 
the  presence  of  magnetic  field  are  usually  performed  at  room 
temperature  using  the  same  techniques  as  those  used  in  the 
absence  of  magnetic  fields.  The  magnetic  field  is  generated  by 
either  electromagnets  or  permanent  magnets.  Richler  et  al.  [44] 
investigated  the  effect  of  magnetically  driven  structure  formation 
on  heat  flux  in  MNFs  on  the  basis  of  thermal  conductivity 
measurements  in  variation  of  an  external  magnetic  field.  For  this 
purpose,  they  developed  an  improved  measuring  device  based  on 
the  plane  heat  source  instead  of  the  standard  hot  wire  method  is 
used  to  enable  both  parallel  and  perpendicular  orientations  of 
magnetic  field  and  heat  flux.  Thermal  conductivity  measurements 
were  carried  out  in  variation  of  strength  and  direction  of  an 
external  magnetic  field  relative  to  heat  flux.  It  was  concluded  that 
unlike  former  experimental  investigations,  for  the  first  time  the 
results  showed  qualitative  consistency  with  theoretical  predic¬ 
tions  for  both  orientations. 

Li  et  al.  [33]  reported  on  the  influence  of  magnetic  field 
strength  and  direction  on  the  transport  properties  of  MNF.  They 
found  that  an  external  magnetic  field  had  remarkable  effects  on 
the  both  thermal  conductivity  and  viscosity  of  the  MNF.  Little 
change  in  the  thermal  conductivity  of  the  MNF  was  found  in  the 
magnetic  field  perpendicular  to  the  temperature  gradient,  irre¬ 
spective  of  the  applied  magnetic  field  strength  and  the  volume 
fraction  of  particles.  The  thermal  conductivity  of  the  MNF 
increased  with  the  strength  of  the  applied  magnetic  field  being 
parallel  to  the  temperature  gradient.  The  reason  for  this  enhance¬ 
ment  was  the  change  of  microstructures  induced  by  the  external 
magnetic  field  in  the  MNF.  They  further  explained  that  when  the 
magnetic  field  was  parallel  to  the  temperature  gradient,  the 
formed  particle  chains  provided  the  more  effectively  bridges  for 
energy  transport  inside  the  MNF  along  the  direction  of  tempera¬ 
ture  gradient  and  enhanced  the  thermal  process  in  the  MNF. 
It  was  revealed  that  the  viscosity  firstly  increased  with  the 
magnetic  field  and  finally  approached  a  constant  as  the  magne¬ 
tization  of  the  magnetic  fluid  achieved  a  saturation  state.  This 
behavior  was  also  explained  by  the  presence  of  chain-like 
structures  in  MNF  with  respect  to  the  influence  of  magnetic  field. 
Philip  et  al.  [22,37]  reported  a  dramatic  enhancement  of  the 
thermal  conductivity  up  to  K/Kf=  4.0  (300%)  in  magnetite-based 
nanofluid,  under  the  influence  of  an  applied  magnetic  field  along 
the  direction  of  heat  flow.  The  reported  thermal  conductivity 
enhancement  was  within  the  predicted  value  for  parallel  mode 
conduction.  The  decrease  in  K  observed  after  the  critical  value  of 
magnetic  field  was  explained  by  the  “zippering”  of  chains.  They 
argued  that  the  cluster  morphology  and  distribution  both  could 
have  impacts  on  thermal  conductivity  enhancement  and  support 
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the  transport  of  heat  through  the  aggregates.  In  their  report,  a 
magnetically  polarizable  nanofluid  (magnetic  nanofluid)  could  be 
used  as  a  reversible  switchable  thermal  fluid,  i.e.,  insulation  to 
high  thermal  conducting  fluid  and  vice  versa,  by  changing  the 
magnetic  field.  Wright  et  al.  [45]  reported  an  innovative  concept 
of  integrating  magnetically  sensitive  metal  or  metal  oxides  in  a 
fluid  with  carbon-nanotubes  in  order  to  increase  the  thermal 
conductivity  of  the  fluid.  It  was  shown  that  the  MNF  with  Ni 
coated  single  wall  nanotubes  was  significantly  enhanced  in  the 
presence  of  magnetic  field.  The  thermal  conductivity  enhance¬ 
ment  observed  was  attributed  to  the  influence  of  magnetic  field 
on  the  microstructure  of  the  fluid.  It  was  indicated  that  under  the 
applied  magnetic  field,  the  Ni  coated  nanotubes  could  form 
aligned  chains  that  favored  to  connect  the  nanotubes,  resulting 
in  the  improved  contacts.  Wensel  et  al.  [46]  investigated  the 
thermal  conductivity  of  nanofluids  with  metal  oxides  nanoparti¬ 
cles  (Fe304  and  MgO)  and  carbon  nanotubes  in  the  presence  of 
magnetic  field.  They  observed  that  a  maximum  thermal  conduc¬ 
tivity  could  be  increased  by  up  to  10%  at  a  rather  low  particle 
weight  fraction  of  0.02  wt%.  The  possible  explanation  to  these 
interesting  results  was  the  aggregation  of  metal  oxide  particles  on 
the  surface  of  nanotubes  by  electrostatic  attraction  and  the 
formation  of  the  chain-structure  along  the  nanotubes.  It  was 
suggested  that  the  investigated  nanofluids  could  be  used  for 
coolant  applications  since  their  viscosity  was  similar  to  that  of 
water.  Nkurikiyimfura  et  al.  [36]  investigated  the  influences  of 
particle  size  and  particle  volume  fraction  on  the  thermal  con¬ 
ductivity  enhancement  of  an  engine  oil  based  magnetite  magnetic 
nanofluid  in  a  magnetic  field.  It  was  shown  that  the  thermal 
conductivity  could  be  enhanced  with  the  smaller  magnetite 
particles  used  in  the  nanofluids  at  a  greater  particle  volume 
fraction.  In  addition,  the  thermal  conductivity  enhancement  was 
found  to  be  related  to  magnetic  field  parallel  to  temperature 
gradient.  They  analyzed  that  the  increase  of  thermal  conductivity 
ratio  was  due  to  the  chain-like  structures  formed  in  MNF  under 
the  influence  of  magnetic  field.  The  thermal  conductivity 
enhancement  observed  for  the  smaller  particles  was  explained 
by  the  super-paramagnetic  behavior  of  the  smaller  particles  and 
their  effect  on  the  aggregates  formation.  Though  the  magnetic 
field  parallel  to  temperature  gradient  exhibited  a  positive  effect 
on  the  thermal  conductivity  of  the  MNF,  a  high  magnetic  field  had 
a  negative  effect  on  other  thermophysical  properties  of  the  fluid 
such  as  rheological  properties  and  thermal  conductivity  itself 
[22,36,37].  This  is  due,  especially,  to  the  undesirable  clumping 
and  zippering  of  chains  under  the  influence  of  high  magnetic 
fields.  This  restriction  calls  for  the  simultaneous  treatment  of 
thermal  conductivity  enhancement  with  other  transport  proper¬ 
ties,  especially,  the  viscosity  (>;).  When  the  thermal  conductivity 
and  magnetoviscous  effects  were  considered,  the  enhancement 
could  be  characterized  by  a  ratio  of  the  thermal  conductivity  to 
the  viscosity.  Philip  et  al.  [47]  addressed  this  issue  and  suggested 
the  tuning  of  the  K/ (/-ratio  rather  than  the  thermal  conductivity 
alone.  Their  report  devoted  to  the  determination  and  character¬ 
ization  of  the  viscosity  and  thermal  conductivity  enhancement  in 
a  stable  magnetic  nanofluid  with  the  particles  size  less  than 
10  nm  as  a  function  of  volume  fraction,  shear  rate,  magnetic  field 
and  time.  Their  results  revealed  that,  without  any  magnetic  field, 
the  increase  of  MNF  viscosity  was  much  greater  than  the  thermal 
conductivity  enhancement.  They  also  noted  that  the  K/^-ratio  of 
MNF  at  0.078  vol%  could  be  tunable  from  0.725  to  2.35  by 
controlling  the  applied  field.  The  MNF  with  tunable  K/i]- ratio 
could  be  used  as  multifunctional  “smart  materials"  for  miniature 
cooling  with  damping  applications.  Gavili  et  al.  [48]  conducted 
experimental  investigation  on  the  thermal  conductivity  of  water 
based  MNF  under  magnetic  field  created  by  the  Helmholtz  coils. 
The  magnetic  field  strength  was  controlled  by  an  electric  current. 


The  saturation  time  and  the  reversibility  of  thermal  conductivity 
were  also  examined  after  the  magnetic  field  was  turned  off.  The 
obtained  results  showed  that  the  MNFs  with  5.0%  volume  fraction 
of  nanoparticles  with  an  average  diameter  of  10  nm  enhanced  the 
thermal  conductivity  more  than  200%  at  maximum  value. 
Furthermore,  the  variation  of  thermal  conductivity  with  tempera¬ 
ture  was  also  shown. 

3.2.  Mechanisms  of  thermal  conductivity  enhancement  and 
theoretical  models 

Though  the  thermal  conductivity  of  MNF  has  been  a  subject  of 
some  investigations,  the  mechanisms  for  interpreting  experimen¬ 
tal  data  both  in  the  absence  and  presence  of  applied  magnetic 
field  are  still  needed.  Among  the  mechanisms  proposed  for 
conventional  nanofluids,  the  Brownian  motion  and  particle  clus¬ 
tering/structure  are  the  two  much-debated  mechanisms  [37,49]. 

3.2.1.  Brownian  motion 

The  Brownian  motion  refers  to  the  seemingly  random  move¬ 
ment  of  particles  suspended  in  liquid  or  gas,  and  the  motion  is 
due  to  collision  with  base  fluid  molecules,  which  makes  particles 
undergo  a  random  walk  motion  [3,50,51].  The  Brownian  motion 
could  contribute  to  the  thermal  conductivity  enhancement  in  two 
ways,  namely,  the  direct  contribution  due  to  motion  of  nanopar¬ 
ticles  that  transports  heat  (diffusion  of  nanoparticles)  and  the 
indirect  contribution  due  to  the  so  called  micro-convection  of  fluid 
surrounding  individual  nanoparticles  [5,51,52].  It  is  hypothesized 
that  the  later  effect  could  set  up  a  current  of  heat  transfer 
between  nanoparticles  and  the  carrier  fluid,  hence  enhancing 
the  thermal  conductivity  of  the  nanofluid.  Recent  studies  rejected 
this  hypothesis,  suggesting  that  the  thermal  conductivity 
enhancement  of  MNF  could  be  explained  by  the  particle  cluster¬ 
ing  [22,37,53].  Philip  et  al.  [22]  showed  that  the  microconvection 
of  the  fluid  medium  around  randomly  moving  nanoparticles  did 
not  affect  the  thermal  conductivity  of  a  nanofluid  and  the 
microconvection  model  overestimated  the  thermal  conductivity 
values.  They  analyzed  the  modest  thermal  conductivity  enhance¬ 
ment  in  the  well-dispersed  particles.  The  diffusion  of  magnetic 
nanoparticles  plays  an  important  role  at  a  low  volume  fraction 
[q>  <2%),  which  could  be  explained  by  the  effective  medium 
(Maxwell)  theory  rather  than  the  effects  associated  with  the 
Brownian  motion  induced  hydrodynamics.  Their  explanation  to 
the  conductivity  enhancement  in  MNF  at  a  high  volume 
fraction  (<p  >  2%)  was  due  to  the  presence  of  dimmers  or  trimmers 
in  the  fluid.  These  results  were  in  a  reasonable  agreement 
with  the  Maxwell-Gannet  model,  especially  at  higher  volume 
fractions. 

Tsai  et  al.  [54,55]  reported  the  effect  of  viscosity  on  the 
thermal  conductivity  of  magnetite  based  nanofluids.  They  noted 
that  the  measured  thermal  conductivity  of  nanofluids  gradually 
approached  a  value  predicted  by  the  Maxwell  equation  with 
increasing  the  viscosity.  Their  results  provided  an  evidence  of 
the  viscosity  of  MNF  did  affect  the  thermal  conductivity.  The 
diffusion  of  magnetic  nanoparticles  could  be  an  important  factor 
that  enhanced  the  thermal  conductivity.  It  was  suggested  that  the 
thermal  conductivity  of  MNF  in  the  absence  of  magnetic  fields 
could  be  predicted  by  combined  models  [56-60]  where  the 
conduction  part  could  be  obtained  from  the  prediction  by  the 
Maxwell  equations 

3.2.2.  Nanoparticle  clustering 

Aggregation  of  particles  into  sparse  clusters,  or  ideally  into 
linear  chains  and  its  influence  on  the  nanofluids  properties  have 
received  a  particular  interest  recently,  and  it  is  expected  to  be  a 
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main  mechanism  of  thermal  conductivity  of  the  nanofluids 
[1,2,51,60-64].  The  presence  of  particle  clusters  in  a  nanofluid 
leads  to  extended  and  highly  conductive  paths  for  the  heat  flow 
and  thus  to  a  quick  heat  transport  along  long  distances  since  heat 
can  be  conducted  much  faster  by  solid  particles  if  compared  to  a 
carrier  fluid  [60].  Bishop  et  al.  [65]  reviewed  different  interactions 
occurring  in  magnetic  nanoparticles  self-assembly  as  well  as  their 
scales.  It  was  shown  that  even  without  the  effect  of  external 
magnetic  field;  the  magnetic-interaction-induced  self-assembled 
aggregation  could  occur  in  MNF.  The  aggregation  structure  of 
magnetic  particles  controls  the  heat  conduction  in  the  MNF.  In  the 
absence  of  an  external  magnetic  field,  the  distribution  of  particles 
in  the  MNF  is  disordered  and  the  thermal  conductivity  of  MNF  is 
isotropic.  When  magnetic  nanoparticles  self-assemble,  they  tend 
to  align  their  magnetic  moments  in  the  direction  of  the  local 
magnetic  field  due  to  the  neighboring  particles  or  the  applied 
fields  [21,66,67].  This  gives  rise  to  the  anisotropy  of  interaction 
enabling  magnetic  nanoparticles  to  form  micron-sized,  one¬ 
dimensional  chains/wires,  rings,  two-dimensional  aggregates  or 
even  three-dimensional  super-lattices.  The  understanding  of  the 
aggregate  formation  mechanisms,  their  distribution/distribution 
morphology  and  their  manipulation  with  an  external  magnetic 
field  are  prerequisite  to  understand  the  mechanisms  of  thermal 
conduction  in  MNF.  Hong  et  al.  [39]  investigated  the  effect  of 
magnetic  nanoparticles  clustering  for  an  ethylene  glycol  based  Fe 
magnetic  nanofluids  without  magnetic  fields.  The  thermal  con¬ 
ductivity  was  determined  as  a  function  of  the  duration  of  the 
application  of  the  ultrasonic  vibration,  which  varied  from  0  min 
(i.e.  no  vibration  applied)  to  70  min.  It  was  indicated  that  thermal 
conductivity  ratio  increased  with  increasing  vibration  time,  and 
the  rate  of  this  increase  became  smaller  for  longer  vibration  time. 
Furthermore,  the  variation  of  thermal  conductivity  with  time 
after  the  application  of  vibration  was  investigated,  and  it  was 
found  that  thermal  conductivity  decreased  with  time.  Variation  of 
average  size  of  clusters  was  also  determined  as  a  function  of  time 
after  the  application  of  vibration,  and  the  cluster  size  increased 
with  time.  As  a  result,  the  size  of  the  clusters  formed  by  the 
nanoparticles  had  a  major  influence  on  the  thermal  conductivity. 
The  effect  of  particle  clustering  on  the  thermal  conductivity  was 
also  claimed  by  Zhu  et  al.  [41  ].  They  found  that  the  clustering  and 
nanoparticle  alignment  was  mainly  responsible  for  the  anoma¬ 
lous  thermal  conductivity  enhancement.  Jiang  et  al.  [68]  prepared 
a  MNF  by  one  step  phase  transfer,  and  investigated  the  aggrega¬ 
tion  effect  on  thermal  conductivity  enhancement  of  the  prepared 
MNF.  The  microstructure  of  the  MNF  was  analyzed  by  a  dynamic 
light  scattering  method,  which  unveiled  the  variation  of  aggre¬ 
gated  configurations  with  particle  concentration  and  time.  The 
observed  wave-like  structure  of  the  measured  thermal  conduc¬ 
tivity  as  a  function  of  particle  concentration  was  explained  as  the 
manifestation  of  the  aggregation-structure  variation.  The  analysis 
for  the  possible  mechanism  of  aggregate  formation  in  MNF  in  zero 
magnetic  fields  revealed  that  the  particle  coating  layers  were  also 
of  importance  in  the  cluster  formation.  They  concluded  that  the 
thermal  conductivity  of  MNF  could  be  manipulated  for  some 
nanofluids  by  changing  a  stabilizer  used  and  thus  controlling  the 
structure  of  the  aggregates.  Philip  et  al.  [22]  performed  an 
experiment  to  prove  the  enhanced  thermal  conduction  through 
percolating  structures  in  MNF.  They  stated  that  the  maximum 
enhancement  was  observed  when  the  chains  like  aggregates 
were  well-dispersed  without  clumping.  The  MNF  could  be  used 
as  a  reversible  switchable  thermal  fluid,  i.e.,  insulation  to  high 
thermal  conductivity  and  vice  versa,  via  the  change  of  the 
magnetic  field. 

Some  attentions  have  been  paid  to  theoretical  and  mathema¬ 
tical  modeling  of  the  thermal  conductivity  enhancement  of  MNF 
both  in  the  absence  and  presence  of  magnetic  fields.  Table  1 


summarizes  the  commonly  used  models  to  validate  experimental 
results. 

Fang  et  al.  [69]  attempted  to  model  the  thermal  conductivity 
of  MNF  in  the  presence  of  magnetic  field.  Their  model  consisted  of 
the  modified  Maxwell-Gannet  model  by  the  introduction  of  an 
anisotropic  constant  that  accounted  for  the  microstructure  of  the 
fluid  under  an  applied  magnetic  field. 

The  obtained  model  was  given  as 

Kx  Kp  +  2Kf  +  2cpx(Kp-I<f) 

Kf  I<P  +  2Kf-(px(Kp-Kf) 


With  <px=(l  +  Cx)<pwhere  Cx  is  the  anisotropic  structure  para¬ 
meter  in  the  x  direction  expressed  as 


i  —  1  j  >  i 


(2) 


(3) 


where  subscripts  i  and  j  are  particle  numbers,  d  is  the  diameter  of 
the  particle,  r,j= r,  —  rtj  =  |  r  ,j|,  e  ,j  =  r  j/r,j,  r,  is  the  position  of 
the  ith  particle,  n  x  is  the  unit  vector  along  the  x  direction.  From 
Eqs.  (2)  and  (3),  the  higher  particle  volume  fraction,  the  shorter 
distance  between  two  particles  r,j-  =  |  r  ,y|,  the  greater  value  of  the 
anisotropic  structure  parameter  Cx  could  be  obtained.  When  the 
chain  like  structure  along  the  x  direction  becomes  more  obvious, 
the  value  of  the  anisotropic  structure  parameter  Cx  becomes 
greater.  The  anisotropic  structure  parameter  Cx  well  characterizes 
the  distribution  of  particles  along  the  x  direction.  Note  that  the 
same  method  can  be  used  to  calculate  the  thermal  conductivity  of 
MNF  along  the  y  and  z  directions.  The  anisotropic  thermal 
conductivity  of  MNF  calculated  by  the  modified  Maxwell- 
Gannet  formula  is  in  agreement  with  those  computed  by  the 
numerical  methods. 

Fu  et  al.  [70]  conducted  theoretical  investigation  on  the 
effective  thermal  conductivity  tensor  for  MNFs  using  a  two-step 
homogenization  method.  The  differential  effective  medium  the¬ 
ory  was  used  to  determine  the  equivalent  thermal  conductivity  of 
magnetizable  nanoparticle  chains  and  then,  the  theory  was 
generalized  to  self-consistent  anisotropic  effective  medium  the¬ 
ory  to  investigate  the  effective  thermal  conductivity  tensors  of 
magnetic  nanofluids.  The  numerical  results  showed  that  the 
aspect  ratio  of  chain-like  aggregated  clusters  plays  an  important 
role  in  enhancement  of  anisotropic  thermal  conductivity.  In 
addition,  the  theoretical  results  on  the  elements  of  thermal 
conductivity  parallel  to  the  fields  and  perpendicular  to  the  fields 
were  in  good  agreement  with  the  experimental  data.  Their  results 
further  confirmed  the  non-monotonic  dependence  of  effective 
thermal  conductivity  on  magnetic  field  strength,  in  accordance 
with  experimental  reports. 

Nkurikiyimfura  et  al.  [71]  investigated  the  effect  of  chain-like 
magnetic  nanoparticles  aggregates  on  the  thermal  conductivity  of 
magnetic  nanofluids  (MNFs)  in  magnetic  field  via  a  proposed 
model  of  thermal  conductivity  of  MNFs.  The  anisotropic  feature  of 
thermal  conductivity  ratio  was  predicted  based  on  the  field 
induced  chain-like  magnetic  nanoparticles  aggregates  in  MNFs. 
It  was  indicated  that  for  a  magnetic  field  parallel  to  temperature 
gradient,  the  thermal  conductivity  component  along  the  magnetic 
field  was  significantly  enhanced  due  to  the  magnetic  field  induced 
chain-like  magnetic  nanoparticles  aggregates.  The  experimental 
data  of  thermal  conductivity  enhancement  along  the  magnetic 
field  direction  were  similar  to  those  predicted  by  the  proposed 
model. 
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4.  Thermomagnetic  convection 


Thermomagnetic  convection  refers  to  a  convective  heat  trans¬ 
fer  that  makes  use  of  the  spatial  gradient  in  the  magnetic 
susceptibility  of  MNF  that  is  produced  in  the  presence  of  a 
temperature  gradient  [72-74].  When  a  MNF  is  exposed  to  a 
non-uniform  magnetic  field  in  the  presence  of  a  temperature 
field,  besides  the  conventional  gravitational  body  force,  the 
varying  susceptibilities  result  in  a  non-uniform  magnetic  body 

force  on  the  MNF  ( Kelvin  body  force),  given  as  f  m  =  / i0(M  ■  V )  B  , 
where  p0  stands  for  the  magnetic  permeability  of  the  vacuum, 

and  M  the  magnetization  and  B  the  magnetic  induction  [75,76], 
When  the  approximation  of  equilibrium  magnetization  with 
magnetic  susceptibility  is  dependent  solely  on  the  local  value  of 
applied  magnetic  field  and  the  density,  magnetic  and  gravita¬ 
tional  forces  both  are  potential.  The  Kelvin  body  force  creates  a 
static  pressure  field  in  the  flow  that  is  symmetric  about  the 
applied  magnetic  field  producing  an  irrotational  force  field.  Such  a 
symmetric  field  does  not  alter  the  velocity  profile  and  the 
convection  inside  the  fluid  cannot  arise  [9[.  The  non-potentiality 
of  the  bulk  forces  (gravitational  and  magnetic  forces)  appears 
only  if  a  fluid  possesses  a  spatial  non-uniformity  of  the  density  p 

and  magnetization  M  due  to  their  dependences  on  temperature 
or  on  the  particle  concentration  cp  [77],  The  later  condition  is 
satisfied  for  non-isothermal  systems  with  an  asymmetric  (or 
biased)  temperature  distribution  about  the  imposed  magnetic 
field.  The  resultant  Kelvin  body  force,  which  is  also  asymmetric  in 
this  case,  creates  a  field  force,  which  leads  to  the  self-organized 
adjective  motion  of  the  MNF  across  isotherms.  The  condition  for 
the  free  convection  to  develop  is  calculated  as  follows  [77]: 


/  =  VTx 


PrPot  ±p0PmH^VH2 


*0 


(4) 


where  p  and  ft T  denote  the  density  of  magnetic  liquid  and  its 
thermal  expansion  coefficient,  respectively.  In  high  magnetic 
fields,  when  magnetization  of  liquid  is  close  to  saturation,  the 
pyromagnetic  coefficient  can  be  simplified  as 

P  =  M(PT  +  Pm)  (7) 

where  /Jm=(l/m)(dm/dr)  denotes  the  relative  pyromagnetic  coef¬ 
ficient  of  dispersed  material.  Under  these  conditions,  the  value  of 
pyromagnetic  coefficient  increases  with  increasing  the  sum 
P=(\dm/dT\  +Pt)  in  the  required  temperature  range  for  the  same 
MNF.  This  condition  suggests  that  the  choice  of  magnetic  material 
should  depend  on  the  thermal  expansion  of  the  carrier  fluid.  For 
liquids  with  small  thermal  expansion,  this  coefficient  of  (dm/dT)  is 
more  important.  In  low  magnetic  fields,  the  linear  law  of  magne¬ 
tization  can  be  assumed  and  the  value  of  pyromagnetic  coefficient 
is  expressed  as 

P  =  M(^+pT  +  2pn^  (8) 

It  is  clear  that  under  low  magnetic  fields,  the  share  of  dispersed 
material  in  pyromagnetic  coefficient  becomes  more  significant  [78]. 
Recent  efforts  have  been  made  to  synthesize  MNF  suitable  for  heat 
transfer  application.  A  special  aspect  is  to  prepare  magnetic  nano- 
particles  with  a  greater  pyromagnetic  coefficient,  i.e.,  at  the  Curie 
temperature  being  close  to  the  operating  range  (considering  the 
boiling  point  of  the  carrier  and  the  stability  of  the  fluid  at  elevated 
temperature).  The  most  used  magnetic  materials  are  the  Zn- 
substituted  spinel-type  ferrites  because  Zn  substitution  alters  their 
magnetic  parameters  in  a  wide  range  of  values  [79-81  ].  This  allows  to 
obtain  ferrite  particles  with  different  thermomagnetic  coefficients. 
Mn0.5Zn0  5Fe204  based  magnetic  nanofluids  are  widely  used  amongst 
due  to  their  low  Neel  temperature  (100  °C),  and  thus  higher  thermo¬ 
magnetic  coefficients  (see  Table  2). 


where  /  denotes  the  vector  sum  of  the  Kelvin  body  force  and  the 
gravitational  body  force;  the  sign  +  +  denotes  the  parallel  and 
anti-parallel  orientation  of  the  magnetic  field  gradient  with 
respect  to  the  gravitational  force,  pm  is  relative  pyromagnetic 
coefficient  of  the  fluid,  and  pT  is  the  relative  volumetric  expansion 
coefficient. 

The  heat  transfer  intensity  is  determined  by  the  Rayleigh 
number  Ra,  which  consists  of  thermomagnetic  and  thermogravi- 
tational  parts: 

Ra  =  Rdr+Ram  =  (^PtPS+ PoPm^l-^)  (5) 


According  to  the  numerical  analysis,  even  in  lab-scale  experi¬ 
ments,  thermal  ( RaT )  and  magnetic  ( Ram )  Rayleigh  numbers  both 
are  greater;  they  can  exceed  the  values  of  105-106,  showing  that 
the  problem  of  convective  heat  transfer  may  be  theoretically 
considered  in  the  framework  of  a  boundary  layer  approximation 
[77].  From  the  expression  of  the  Rayleigh  number,  it  is  clear  that 
the  intensity  of  thermomagnetic  convection  and  therefore  the 
efficiency  of  devices  using  MNF  as  heat  transfer  media  is  provided 
not  only  by  the  magnetic  field  and  the  temperature  distributions 
fields  but  also  by  the  properties  of  MNF  with  the  thermophysical 
properties  and  the  pyromagnetic  coefficient,  i.e.,  the  degree  of 
dependence  of  magnetization  on  temperature.  In  general,  the 
pyromagnetic  coefficient  of  a  MNF  is  given  by  [78] 


P  =  PTM+ 


/  8M\ 

dm 

5M 

\dm) 

df 
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with 
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1  dp 
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4.1.  Experimental  investigations  on  thermomagnetic  convection 

Shuchia  et  al.  [82]  investigated  the  effect  of  magnetic  field 
position  on  heat  transfer  and  driving  force  characteristics  of  a 
mixture  of  MNF  with  an  organic  fluid  introduced  in  a  heat 
transport  device.  Their  results  indicated  that  the  heat  transfer 
capability  and  the  magnetic  driving  force  were  improved  by 
magnetic  fields  when  the  magnetic  field  was  applied  on  the 
entrance  of  the  heated  region.  Blums  et  al.  [77]  investigated  the 
heat  transfer  from  a  nonmagnetic  cylinder  to  a  MnxZn(1_x)Fe204 
tetradecane  temperature  sensitive  magnetic  nanofluid,  under  the 
influence  of  transverse  laminar  free  convection  in  the  presence  of 
uniform  and  non-uniform  magnetic  fields  directed  transversally  to 
the  axis  of  the  cylinder.  Their  results  confirmed  the  additive  action 
of  thermo-gravitation  and  thermomagnetic  body  forces  on  the  heat 
transfer  intensity  in  ferrofluids.  The  heat  transfer  enhancement 
was  attributed  to  the  properties  of  MNF  and  the  magnetic  field 
gradient.  They  explained  that  it  was  possible  to  use  temperature 
sensitive  magnetic  nanofluids  to  achieve  a  magnetic  Rayleigh 
number  ( Ram )  that  exceeds  the  thermo-gravitational  Rayleigh 
number  (Rar).  They  also  found  that  a  greater  magnetic  Rayleigh 
number  and  a  higher  heat  transfer  could  be  achieved  when  the 
magnetic  field  gradient  was  increased.  Zablockis  et  al.  [83]  con¬ 
ducted  experimental  and  numerical  investigations  of  thermomag¬ 
netic  convection  in  a  heated  cylinder  under  the  magnetic  field  in  a 
solenoid.  They  found  that  the  heat  transfer  under  a  magnetic  field 
could  be  increased.  However,  the  solenoid  used  showed  to  be 
impractical,  especially  for  generating  sufficiently  strong  magnetic 
fields.  Yamaguchi  et  al.  [84]  conducted  an  investigation  on  a  Mn- 
Zn  ferrite  alkyl-naphthalene  based  magnetic  nanofluid  in  a  cubic 
cavity.  The  magnetic  field  and  magnetization  were  affected  by 
temperature,  and  the  greater  magnetization  and  magnetic  forces 


I.  Nkurikiyimfura  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  21  (2013)  548-561 


555 


Table  2 

Summary  of  theoretical  and  numerical  investigations  on  thermomagnetic  convection. 


Author 

Geometry  and 
Dimension 

Magnetic  field 

MNF 

Model 

Remarks 

Zablockis 

Cylindrical 

Non-uniform  constant  magnetic 

Temperature 

Lattice 

The  HTC  increased  up  to  sevenfold  within  the  studied  range 

et  al.  [89] 

domain 

field  of  a  Solenoid  placed  in  a 

sensitive  MNF 

Boltzmann 

of  parameters 

dimensions 

2x6  cm 

hollow  area  inside  the  domain 

model 

Yamaguchi 

Square  cavity 

Uniform  magnetic  field  generated 

Alkyl-naphthalene 

Lattice 

The  HTC  increased  with  the  magnetic  field  strength 

et  al.  [91] 

by  vertical  electromagnetic 

based  Mn-Zn 

Boltzmann 

ferrite  MNF 

model 

Yamaguchi 

Partitioned 

Alkyl-naphthalene 

Single  model 

Heat  transfer  enhancement  is  affected  by  improving  heat 

et  al. 

rectangular  box 

based  Mn-Zn 

transfer  characteristics  at  higher  Rayleigh  number  when  a 

[104] 

ferrite  MNF 

strong  magnetic  field  is  imposed 

Yamaguchi 

Cubic  cavity 

Uniform  magnetic  field  generated 

Mn-Zn  ferrite  alkyl- 

Lattice 

Heat  transfer  characteristic  of  the  used  MNF  was  enhanced 

et  al.  [90] 

by  permanent  magnets 

naphthalene  based 

Boltzmann 

when  the  magnetic  field  was  applied 

MNF 

model 

Mukhopad- 

Two- 

Magnetic  field  created  by  a  line 

— 

Single  phase 

The  Nusselt  number  scales  with  Raj25  and  thermomagnetic 

hyay 

dimensional 

dipole 

model 

heat  transfer  increases  when  the  length  scale  decreases 

et  al.  [82] 

enclosure 

(scaling 

analysis) 

Zablotsky 

Rectangular  cell 

Permanent  magnet  attached  to  the 

Tetradecane 

Single  phase 

Thermomagnetic  convection  exceeds  the  thermo- 

et  al.  [96] 

dimensions: 

cell 

temperature 
sensitive  Mn-Zn 

model 

gravitational  convection: 

100  x  15  mm 
height:  150  mm 

based  MNF 

and  a  higher  enhancement  is  obtained  when  the  magnet  is 
placed  near  the  warmer  source 

Jafar  et  al. 

Cylinder 

_ 

Kerosene  based 

Mixture 

The  heat  transfer  enhancement  is  function  of  magnetic 

[105] 

MNF 

(two-phase) 

strength  and  direction.  The  heat  transfer  enhancement 

Dimensions: 

3.5  x  75  mm 

model 

decreases  with  the  increase  of  particle  size 

Jue  [106] 

2-D  square  cavity 

Single  phase 

The  magnetic  nanofluid  flow  and  heat  transfer  are  adjusted 

model 

not  only  by  the  magnet  strength  but  also  by  a  proper  choice 

Dimensions: 
w  =  0.05, 

of  its  position 

a  =  3.475  mm, 
b  =  3.475  mm, 

L  =  5.5  mm 

Strek  and 

Channel  between  Dipole  located  below  the  channel  — 

Single  model 

External  magnetic  field  has  influence  on  both  on  flow  and 

Joseph 

two  parallel 

(COSMOL 

heat  transfer  and  the  convection  process  is  similar  to  the 

[107] 

plates 

Dimensions: 

h  —  0.02  m ; 

L  =  0.2m 

codes) 

thermo-  gravitational  convection 

Jafar  et  al 

Cylindrical 

Mixture 

The  magnetic  soret  effect  can  be  higher  than  the 

[102] 

geometry 

model 

conventional  one  and  its  strength  depends  on  the  magnetic 
field  magnitude 

Xuan  et  al. 
[108] 

Micro-channel  — 

Lattice 

Boltzmann 

model 

The  heat  transfer  enhancement  is  related  to  the  orientation 
magnetic  field  and  a  higher  enhancement  is  obtained  when 
the  magnetic  field  gradient  is  parallel  with  the  incoming 
stream. 

were  observed  in  the  regions  near  the  upper  wall  and  center  inside 
the  cavity  rather  than  in  the  region  near  the  bottom  and  side  walls. 
The  analysis  of  the  fluid  flow  showed  a  weak  flow  roll  inside  cavity 
under  the  magnetic  force,  which  brought  the  lower  temperature 
fluid  downward  in  the  center  region,  and  streamed  the  higher 
temperature  fluid  upward  along  the  regions  near  the  side  walls. 
They  also  found  that,  with  the  magnetic  field  imposed,  the  heat 
transfer  inside  the  cavity  was  increased,  compared  to  that  without 
the  magnetic  field,  and  the  heat  transfer  increased  further  with 
increasing  the  strength  of  the  magnetic  field.  Yamaguchi  et  al.  [85] 
also  performed  the  numerical  and  experimental  investigations  on 
heat  transfer  characteristics  of  a  temperature  sensitive  magnetic 
nanofluid  filled  in  a  cubic  container  with  a  heat  generating  square 
cylinder  stick  inside.  Their  results  showed  that,  regardless  of  the 
heat  generating  object  sizes,  the  heat  transfer  characteristic  of  the 


MNF  was  enhanced  when  the  magnetic  field  was  applied.  How¬ 
ever,  the  heat  transfer  process  became  poor  as  the  size  of  the  inside 
heat  generating  object  increased.  This  poor  heat  transfer  process 
was  due  to  the  decrease  of  the  space  where  the  fluid  flowed  with 
the  increase  of  the  heat  generating  object  size.  The  experimental 
data  were  predicted  by  the  numerical  result,  disclosing  the  more 
flow  details  of  the  natural  convection  of  the  MNF  inside  cavity. 

Some  previous  work  focused  on  the  synergy  between  the 
magnetic  field  and  temperature  gradient  and  thermomagnetic 
convection  was  analyzed  irrespective  to  the  effect  of  magnetic 
field  on  the  thermophysical  properties  of  the  used  MNF.  However, 
recent  studies  have  shown  that  the  thermomagnetic  convection 
can  be  strongly  affected  by  the  externally  applied  magnetic  field 
through  the  thermophysical  properties  of  the  fluid  [86,87].  Engler 
et  al.  [86]  investigated  the  magnetoviscous  effect  on  the 
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thermomagnetic  convection.  Under  certain  circumstances,  the 
magnetoviscous  effect  can  affect  the  thermomagnetic  convection. 
The  increase  of  magnetoviscous  effect  was  attributed  to  the 
presence  of  coarse  particles  in  the  studied  samples.  It  was 
concluded  that  for  the  thermomagnetic  convection,  the  magne¬ 
toviscous  effect  could  be  increased  in  the  MNF,  which  would 
affect  the  experimental  results  of  the  threshold  of  convection.  The 
influence  of  magnetoviscous  effects  on  the  thermomagnetic  con¬ 
vection  was  also  confirmed  in  a  recent  theoretical  investigation 
on  the  onset  of  Bernard-Mangaroni  thermomagnetic  convection 
in  a  horizontal  magnetic  nanofluid  layer  by  Nanjundappa  et  al. 
[88].  Also,  Lajvardi  et  al.  [89]  addressed  the  forced  thermomag¬ 
netic  convection  heat  transfer  of  a  water  based  magnetic  nano¬ 
fluid  with  nanoparticles  of  Fe304  (10-nm)  in  a  heated  copper  tube 
for  laminar  flow.  They  observed  a  significant  enhancement  in  heat 
transfer  coefficient  under  the  influence  of  an  applied  magnetic 
field  and  magnetic  nanoparticle  volume  fraction.  The  observed 
heat  transfer  enhancement  was  attributed  to  the  thermophysical 
properties  of  MNF  such  as  thermal  conductivity  or  heat  capacity 
under  the  applied  magnetic  field. 

4.2.  Mathematical  modeling  and  simulation 

Several  numerical  investigations  have  been  conducted  on 
different  geometries  such  as  the  Rayleigh-Bernard  configuration, 
rectangular  enclosure,  cylinder,  cube,  partitioned  cavity,  etc. 
[90,91].  The  heat  transfer  processes  and  fluid  flow  in  these 
geometries  have  been  simulated  in  the  presence  of  static  mag¬ 
netic  field  and  more  recently  in  modulated  magnetic  fields.  The 
commonly  used  mathematic  models  can  be  summarized  in  two 
categories: 

1.  Single  phase  (homogeneous)  model.  MNF  is  a  multiphase  mix¬ 
ture  of  magnetic  nanoparticles  and  carrier  fluid.  However,  due 
to  the  restricted  magnetic  nanoparticle  size,  the  assumptions 
of  thermal  equilibrium  conditions  and  negligible  motion  slip 
between  particles,  MNF  is  usually  considered  as  a  conventional 
single  phase  fluid  with  average  physical  properties  of  individual 
phase  that  exhibit  magnetic  behavior.  In  this  case,  the  conven¬ 
tional  hydrodynamic  equations  can  be  modified  to  account  for 
the  Kelvin  body  force,  the  so  called  ferro-hydrodynamic  equations, 
and  coupled  with  magnetization  equation  and  Maxwell’s 
equations  used  to  model  and  characterize  the  fluid  flow  and 
heat  transfer  processes.  Mahmoudi  et  al.  [92]  analyzed  the 
natural  convection  for  a  two-dimensional  triangular  enclosure 
with  partially  heated  from  below  and  cold  inclined  wall  filled 
with  nanofluid  in  presence  of  magnetic  field.  The  governing 
equations  were  solved  by  a  finite  volume  method.  The  flow 
pattern,  isotherms  and  average  Nusselt  number  were  presented 
for  0  <  Ha  <  100,  104  <  Ra  <  107,  0  <  q>  <  0.05  and  six  cases  that 
were  made  by  location  of  heat  sources.  The  results  showed  that  in 
presence  of  magnetic  field  the  flow  field  was  suppressed  and  heat 
transfer  decreased.  Furthermore,  it  was  observed  that  a  maximum 
reduction  of  average  Nusselt  number  in  a  high  value  of  Ha 
occurred  when  Ra=  106.  It  was  found  that  the  nanoparticles  were 
more  effective  when  Ra  =104  where  conduction  was  more  pro¬ 
nounced.  Sheikholeslami  et  al.  [93]  investigated  the  natural 
convection  in  a  concentric  annulus  between  a  cold  outer  square 
and  heated  inner  circular  cylinders  in  the  presence  of  static  radial 
magnetic  field  using  the  lattice  Boltzmann  method.  The  inner  and 
outer  cylinders  were  maintained  at  constant  uniform  tempera¬ 
tures,  and  it  was  assumed  that  all  walls  could  insulate  the 
magnetic  field.  The  numerical  investigation  was  carried  out  for 
different  governing  parameters,  i.e.,  the  Hartmann  number,  nano- 
particles  volume  fraction  and  the  Rayleigh  number.  The  effective 
thermal  conductivity  and  viscosity  of  nanofluids  were  calculated 


using  the  Maxwell-Garnetts  (M-G)  and  Brinkman  models,  respec¬ 
tively.  Also,  the  multi-distribution-function  (MDF)  model  was 
used  for  the  uniform  magnetic  field  effects  simulation.  The  results 
revealed  that  the  average  Nusselt  number  was  an  increasing 
function  of  nanoparticle  volume  fraction  as  well  as  the  Rayleigh 
number,  while  it  was  a  decreasing  function  of  the  Hartmann 
number. 

2.  Mixture  model  (two-phase  model).  The  mixture  model  uses  a 
single  fluid  approach,  and  it  is  an  intermediate  between  the 
single  phase  approximation  and  full  set  of  equations  governing 
the  dynamic  of  multiphase  flow  [94].  For  this  reason,  some 
factors  such  as  friction  between  magnetic  nanoparticles  and 
the  carrier  fluid,  the  Brownian  diffusion,  sedimentation  and 
dispersion  are  included  in  the  model.  The  inclusion  of  the  later 
factors  in  the  model  provides  a  possibility  of  understanding 
the  function  of  each  phase  in  the  heat  transfer  process  and  the 
problem  related  to  the  stability  of  MNF  under  both  magnetic 
and  gravitational  forces  [95,96]. 

Mousavi  [97]  used  CFD  techniques  to  investigate  the  effects  of 
single  phase  approximation  and  mixture  model  on  the  flow 
behavior  and  heat  transfer  of  a  kerosene  based  magnetic  nano¬ 
fluid  in  a  cylindrical  geometry.  It  was  found  that  at  a  great  particle 
volume  fraction,  the  mixture  model  approach  is  more  effective 
than  single  phase  approach,  and  the  former  approach  is  conve¬ 
nient  to  study  the  effects  of  aggregation  and  particle  size  on 
hydrodynamics  of  the  system.  Using  the  mixture  model,  they 
obtained  that  with  coarse  magnetic  particles,  the  heat  transfer 
could  decrease  and  the  Rayleigh  rolls  would  not  be  observed. 

Xuan  et  al.  [98]  investigated  the  flow  and  thermal  processes  of 
MNF  flowing  through  a  micro-channel.  By  altering  the  orientation 
and  magnitude  of  an  external  magnetic  field,  they  could  control 
the  flow  and  thermal  processes  of  MNF.  They  also  found  that  the 
effect  of  the  external  magnetic  field  intensity  was  dependent 
upon  its  magnitude,  its  synergistic  relation  with  the  flow  direc¬ 
tion  of  the  main  stream,  the  temperature  gradient,  and  gravity  in 
the  gravitational  field.  They  concluded  that  for  the  purpose  of 
heat  transfer  enhancement,  the  most  remarkable  effect  could  be 
achieved  when  the  orientation  of  the  magnetic  field  gradient  is 
parallel  with  the  incoming  stream.  They  analyzed  that  a  magnetic 
field  gradient  opposite  to  the  main  stream  would  suppress  a  heat 
exchange  between  the  MNF  and  the  solid  surface.  Jafari  et  al.  [99] 
used  the  mixture  model  to  simulate  the  heat  transfer  phenomena 
in  a  kerosene  based  magnetic  nanofluid  flowing  in  a  cylindrical 
geometry.  It  was  indicated  that  the  increase  of  magnetic  nano¬ 
particle  size  could  cause  the  generation  of  colloids  in  the  system, 
which  had  a  detrimental  effect  on  the  heat  transfer  of  the  system. 
Jafari  et  al.  [96]  investigated  the  parameter  effects  on  the 
thermomagnetic  convection  using  the  Taguchi  technique.  The 
thermomagnetic  convection  in  a  cylindrical  geometry  depended 
on  some  parameters  like  the  particle  diameter,  particle  volume 
fraction,  temperature  difference  in  the  magnetic  nanofluid  layer, 
the  magnetic  field  magnitude  parallel  to  temperature  gradient 
and  the  aspect  ratio  of  the  geometry  as  well.  They  reported  that 
the  magnetic  soret  effect  could  be  greater  than  the  conventional 
one  and  its  strength  depended  on  the  magnetic  field  strength. 
Their  results  were  in  agreement  with  the  experimental  data  by 
Volker  et  al.  [100].  Aminfar  et  al.  [101]  used  the  two  phase 
mixture  model  and  the  control  volume  technique  to  investigate 
the  hydro-thermal  characteristics  of  water  based  magnetite 
nanofluid  contained  in  a  vertical  rectangular  duct.  The  duct  was 
exposed  to  a  non-uniform  transverse  magnetic  field  generated  by 
an  electric  current  going  through  a  wire  located  parallelly  under 
the  duct.  The  results  showed  that  the  use  of  the  aforementioned 
magnetic  field  increases  the  Nusselt  number  and  friction  factor 
and  also  creates  a  pair  of  vortices  that  enhances  heat  transfer  and 
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Table  3 

Existing  thermomagnetic  convection  heat  transfer  correlations. 


Reference 

Geometry 

Correlations 

Ashouri  et  al.  [113] 

Square  cavity 

Nu  =  0.03064^-^7— r  ^2278x 

pr0.01S18/(Pr-0.5)(.aia2)703(l-04,)Ka^ 

Mukhopadhyay  [82] 

Rectangular  cavity 

Nu~f 

Moraveji  and  Hejazian  [114] 

Horizontal  tube 

Nu  =  0.00248Re1  03  Pr°  5  ( 1  +tp  )47  5 

Nu:  Nuselt  number,  Nu:  overall  Nuselt  number,  h:  the  heat  transfer  coefficient  (HTC),  U]  and  a2  are  the  constants. 


prevents  sedimentation  of  nano-particles.  Furthermore,  unlike 
the  axial  non-uniform  magnetic  field,  the  increase  of  the  Nusselt 
number  for  the  transverse  magnetic  field  was  considerable  in  all 
lengths  along  the  duct,  and  it  was  also  concluded  that  with 
increasing  the  Reynolds  number,  the  effect  of  the  transverse 
non-uniform  magnetic  field  on  the  Nusselt  number  was  more 
than  that  of  the  axial  non-uniform  magnetic  field.  Suslov  [91] 
investigated  the  linear  stability  of  MNF  between  two  vertical 
differentially  heated  plates  placed  in  a  uniform  external  magnetic 
field  perpendicular  to  the  plates  and  presented  a  complete 
stability  diagram  for  two  and  three-dimensional  disturbances. 
It  was  shown  that  the  thermogravitational  and  magnetic  effects 
were  the  two  mechanisms  of  appearance  of  three  stability  modes. 
A  work  with  the  later  modes  showed  that,  depending  on  the 
governing  parameters,  the  instability  patterns  were  vertical  sta¬ 
tionary  magneto-  convection  rolls  and/or  vertical  or  obliquely 
counter-propagating  thermogravitational  or  thermomagnetic 
waves  [102].  Suslov  et  al.  [102]  used  the  disturbance  energy 
analysis  to  highlight  the  physical  mechanisms  driving  convection 
in  MNF.  They  found  that  the  increase  of  magnetic  effect  was 
manifested  by  the  existence  of  thermomagnetic  waves,  which 
gradually  replaced  thermogravitational  buoyance-driven  waves. 
Krakov  et  al.  [103]  investigated  the  influence  of  porous  media  and 
uniform  magnetic  field  on  the  thermal  convection  in  magnetic 
fluid.  In  porous  square  cavity,  the  competition  between  gravity 
convection  and  thermomagnetic  convection  mechanisms  could 
lead  to  a  complicated  dependence  of  the  heat  flux  through  the 
cavity  in  magnetic  field.  The  increase  of  magnetic  field  could 
enhance  and  depress  heat  transfer.  It  was  also  indicated  that  the 
dependence  of  the  Nuselt  number  on  the  Rayleigh  number  could 
be  complicated,  leading  to  a  hysteresis  in  some  cases.  Ashouri 
et  al.  [104]  analyzed  the  thermomagnetic  convection  heat  transfer 
in  a  two-dimensional  square  cavity  by  a  semi-implicit  finite 
volume  method.  In  their  investigation,  the  side  walls  of  the  cavity 
were  heated  at  different  temperatures,  the  top  and  bottom  walls 
were  isolated,  and  a  permanent  magnet  was  located  near  the 
bottom  wall.  They  noted  that  in  the  absence  of  magnetic  field,  the 
heat  transfer  was  only  affected  by  conduction.  However,  the 
convective  flow  and  heat  transfer  increased  with  increasing 
magnetic  field  intensity  when  a  magnetic  field  was  applied.  Based 
on  the  numerical  analysis,  they  introduced  a  general  correlation 
for  the  overall  Nusselt  number  on  the  side  walls  in  a  wide  range 
of  effective  parameters,  viz.  l<Pr<103,  2  x  106  <Ram  <  10n, 
0  <  Ma  <  1 02,  0  <  t;  <  1 02,  0<Tr<0.3,  0.25  <L*<  1.0,  0<W*< 
0.1,  0.02  <  G*  <  0.1,  where  Pr  is  the  Prandlt  number,  Ram  the 
magnetic  Rayleigh  number,  Ma  the  thermomagnetic  number,  t] 
the  viscosity,  Tr  the  rational  temperature,  G*  the  gap  between  the 
surface  of  permanent  magnet  and  cavity  in  dimensionless  form, 
I*,  W*  denote  the  length  and  width  of  the  permanent  magnet  in 
dimensionless  form  respectively  (see  Table  3).  The  maximum 
error  produced  via  this  correlation  was  approximately  6%.  Mor- 
aveji  et  al.  [105]  presented  the  computational  fluid  dynamics 
(CFD)  method,  with  a  single  phase  approach,  to  determine  the 


effects  of  nanoparticle  concentration  and  flow  rate  on  the  con¬ 
vective  heat  transfer  and  friction  factor  of  water  based  magnetite 
nanofluid  flowing  through  a  plain  copper  tube  in  turbulent 
regime  with  different  Reynolds  numbers  (3000  <  Re  <22000). 
Magnetite  nanoparticles  with  the  average  diameter  of  36  nm, 
suspended  in  water  as  a  base  fluid  with  four  particle  concentra¬ 
tions  of  0.02,  0.1,  0.6  vol%  were  used.  Applying  the  modeling 
results,  they  found  two  relations  to  estimate  the  Nusselt  number 
and  friction  factor,  based  on  the  dimensionless  numbers.  The 
results  showed  that  the  modeling  data  were  in  good  agreement 
with  the  experimental  data.  The  maximum  error  was  around  10%. 

Recent  studies  analyzed  the  thermomagnetic  convection  in  the 
presence  of  modulated  magnetic  field  [106-108].  Kaloni  et  al. 
[107]  performed  a  theoretical  investigation  of  the  convective 
instability  problem  in  the  thin  horizontal  layer  of  a  MNF  heated 
from  below  and  under  alternating  magnetic  fields.  Lange  et  al. 
[106]  considered  the  thermomagnetic  convection  in  a  horizontal 
layer  of  magnetic  field  subjected  to  a  vertical  temperature 
gradient  and  a  spatially  modulated  magnetic  field.  Their  investi¬ 
gations  revealed  that,  in  contrast  to  the  purely  thermal  driven 
system,  the  nonzero  flow  field  of  the  initial  state  was  character¬ 
ized  by  a  2-vortex  structure.  They  also  stated  that  the  later  state 
was  unknown  in  the  classical  Rayleigh-Benard  configuration,  and 
showed  a  potential  of  spatial  modulations  of  the  external  driving 
to  open  new  horizons  in  the  field  of  pattern  formation  with  soft 
magnetic  substances.  Engler  et  al.  [109]  conducted  theoretical 
and  experimental  investigations  on  the  onset  of  convection  under 
the  influence  of  time-modulated  magnetic  field.  The  shift  in  the 
onset  of  convection  depended  on  the  frequency  of  the  external 
magnetic  field.  Matura  et  al.  [108]  investigated  the  influence  of  a 
time-periodic  and  spatially  homogeneous  magnetic  field  on  the 
linear  stability  properties  and  on  the  nonlinear  response  of  a 
ferrofluid  layer  heated  from  below  and  above.  They  observed  that 
the  stability  of  the  conductive  state  and  particularly  the  type  of 
response,  which  can  be  harmonic  or  sub-harmonic,  was  deter¬ 
mined  by  the  system  parameters.  The  stability  boundary  of  the 
conductive  state  in  the  high-frequency  limit  coincided  with  the 
stationary  stability  boundary  when  a  mean  magnetic  Rayleigh 
number  was  used.  The  stability  boundary  for  low-frequency 
modulation  was  shifted  in  a  way  that  the  conductive  state  got 
stabilized.  They  noted  that  sub-harmonic  response  was  not 
typical  for  ferrofluids  because  of  their  high  Prandtl  numbers. 
However,  in  low-Prandtl  number  simulations,  they  found  non¬ 
linear  relaxed  sub-harmonic  convective  states  as  predicted  by  the 
linear  analysis.  They  further  investigated  the  nonlinear  response 
for  heating  from  below  and  above  for  a  large  span  of  modulation 
frequencies.  For  a  high-frequency  modulation,  the  dynamics  was 
nearly  averaged,  and  the  order  parameters  approach  the  values, 
which  have  a  driving  with  the  corresponding  mean  magnetic 
Rayleigh  number.  In  the  case  of  heating  from  above,  a  change  in 
the  sense  of  rotation  of  the  convection  rolls  was  observed. 
They  also  noted  that,  for  a  low-frequency  modulation,  the 
numerical  noise  could  affect  the  pattern  selection.  This  was  the 
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case  when  the  convection  amplitudes  cannot  be  below  the  noise 
level,  depending  on  the  numerical  accuracy,  in  the  under-critical 
phase  of  driving.  The  simulations  showed  that  for  rather  low 
modulation  frequencies  the  oscillation  profiles  approach  the 
stationary  curves. 


5.  Applications 

MNF  can  be  used  as  a  coolant  (thermal  management  applica¬ 
tions)  and/or  a  heat  transfer  medium  in  energy  conversion 
systems.  The  later  application  is  in  the  embryonic  stage.  The  use 
of  MNF  in  thermal  management  is  currently  used  in  some 
commercial  applications  such  as  loud  speakers  cooling  [21].  The 
most  recently  developed  applications  will  be  reviewed  as  follows. 

5.2.  Energy  convention  devices 

Shimada  [110]  analyzed  the  flow  and  heat  transfer  characteristics 
in  a  parallel  duct-type  energy  conversion  device  with  a  MNF  under  a 
non-uniform  magnetic  field  using  the  theoretical  equations  based  on 
the  heat  conduction  theory.  The  results  showed  that  the  heat 
conduction  model  was  capable  of  predicting  the  temperature  and 
pressure  distributions  quantitatively  for  the  performance  of  various 
types  of  magnetic  field.  Lian  et  al.  [Ill]  established  a  mathematical 
model  for  the  prediction  of  flow  and  heat  transport  features  of  the 
temperature-sensitive  magnetic  fluid  and  the  design  of  an  automatic 
energy  transport  device  based  on  the  thermomagnetic  effect.  Their 
results  revealed  that  a  stable  circulation  flow  could  be  maintained  in 
a  loop-shape  channel  in  the  presence  of  a  proper  external  magnetic 
field  and  temperature  gradient  of  the  magnetic  fluid.  In  addition, 
some  studies  on  the  factors  affecting  the  device  performance,  the 
magnetic  field  strength  and  the  fluid  temperature  difference  between 
the  heating  section  and  cooling  section  were  preponderant.  It  was 
indicated  that  the  developed  device  was  not  sufficiently  adequate  to 
be  used  for  heat-to-power  energy  conversion,  and  it  could  be  rather 
used  as  an  automatically-cooling.  The  developed  model  was  validated 
with  the  experimental  data.  The  established  model  could  be  used  to 
simulate  the  performance  of  devices  exploiting  the  thermomagnetic 
effect  for  the  design  and  evaluation  purposes.  Lian  et  al.  [112] 
reported  the  performance  of  automatic  energy  transport  in  cooling 
devices  based  on  the  thermomagnetic  effect  of  a  temperature 
sensitive  magnetic  nanofluid.  The  investigated  devices  consisted  of 
ordered  loop  of  permanent  magnets,  heat  sources,  heat  sinks  and 
temperature  sensitive  magnetic  fluid  assembled  into  an  automatic 
energy  transport  device.  By  adjusting  the  external  magnetic  field 
and/or  temperature  gradient  field  inside  the  magnetic  nanofluid,  they 
could  control  the  energy  transport  process  of  designed  devices.  In 
addition,  the  constructed  device  showed  a  self  regulating  feature  that 
the  flow  velocity  of  the  MNF  increased  with  the  increase  of  the  heat 
load  and  vice  versa.  Their  results  also  showed  that  the  performance  of 
automatic  energy  transport  systems  was  related  to  the  structure  of 
constructed  loop. 

5.2.  Thermomagnetic  convection  based  cooling 

The  thermomagnetic  convection  based  on  cooling  is  one  of  the 
thermal  applications  of  MNF.  As  mentioned  above,  the  application 
of  an  external  magnetic  field  on  a  MNF  with  varying  susceptibil¬ 
ities  gives  a  non-uniform  magnetic  body  force  (i.e.,  Kelvin  body 
force),  which  leads  to  the  thermomagnetic  convection.  Since  most 
of  MNFs  exhibit  a  superparamagnetic  behavior,  they  do  not  obey 
the  Curie-Weiss  law,  thus  becoming  less  magnetic  at  higher 
temperatures.  Recent  development  of  thermomagnetic  cooling 
devices  is  mainly  motivated  by  their  great  potential  application 


for  small  scale  cooling  devices  such  as  in  miniature  micro-scale 
electronic  devices. 

Zablotsky  et  al.  [90]  investigated  a  possibility  of  technical 
application  of  surface  cooling  based  on  thermomagnetic  convec¬ 
tion.  They  found  that  the  cooling  effect  of  75  W/cm2  reached  in 
their  setup  could  meet  only  low  technical  demands.  Their  results 
also  confirmed  that  that  a  significant  augment  of  the  cooling 
intensity  by  thermomagnetic  convection  could  be  obtained  only 
in  a  case  if  the  heat  source  was  located  into  the  region  of  a 
maximal  magnetic  field  intensity,  when  the  efficiency  of  power 
generation  cycle  reached  its  maximum.  A  high  intensification  of 
the  heat  transfer  observed  when  the  magnet  was  placed  in 
proximity  of  the  warmer  source  was  in  agreement  with  the 
previous  studies  by  Rosensweig  [113]  and  Jue  [114].  Li  et  al. 
[115]  reported  the  operating  characteristics  of  a  miniature  auto¬ 
matic  cooling  device  at  different  parameters.  They  found  that  for 
the  design  of  a  cooling  device  with  a  MNF  as  a  coolant,  the 
following  points  should  be  concerned: 

i.  To  select  a  suitable  magnetic  fluid  with  a  greater  pyromagnetic 
coefficient,  a  lower  Curie  point,  a  higher  saturation  magnetiza¬ 
tion,  a  lower  viscosity  and  a  higher  boiling  point; 

ii.  To  exert  a  stronger  intensity  of  the  external  magnetic  field 
and/or  the  field  gradient,  if  possible;  in  order  to  maintain  a 
synergy  between  the  magnetic  field  and  the  temperature 
gradient  field. 

Xuan  et  al.  [116]  reported  the  practical  design  of  a  liquid 
cooling  device  based  on  the  thermomagnetic  effect  for  electronic 
cooling,  in  which  the  waste  heat  generated  from  the  electronic 
elements  was  used  as  a  useful  power  to  drive  the  magnetic  fluid 
flow  and  to  transport  the  waste  heat  to  a  distant  radiator.  They 
found  that  the  thermomagnetic  cooler  with  a  magnet  located 
nearby  the  heat  source  could  realize  a  better  cooling  performance 
and  a  steadier  start-up  procedure.  They  also  emphasized  on  the 
advantages  and  disadvantages  of  thermomagnetic  cooling  for 
electronic  cooling  applications. 

5.3.  Thermal  conduction  based  and  smart  cooling 

A  recent  work  demonstrated  that  the  thermal  conductivity  of 
MNF  could  be  tuned  by  an  externally  applied  magnetic  field  [22]. 
These  findings,  developed  at  the  India  Gandhi  Center  of  Atomic 
Research,  were  used  in  development  of  “smart”  cooling  devices 
that  took  advantage  of  the  large  thermal  conductivity  enhance¬ 
ment  and  the  conduction  to  remove  the  heat  from  a  hot  device 
immersed  in  MNF.  Philip  et  al.  [53]  developed  a  multifunctional 
device  that  used  MNF  with  tunable  thermal  conductivity  to 
viscosity  ratio  to  remove  heat  and  arrested  vibrations  (damper), 
simultaneously.  They  concerned  that  the  development  of  such 
devices  could  have  applications  in  micro-fluidic  devices,  micro- 
and  nano-electromechanical  systems  (MEMS  and  NEMS)  and 
other  miniature  devices.  However,  these  thermal  conductivity 
based  cooling  devices  are  still  in  their  embryonic  stage  and  need  a 
further  development  for  practical  applications. 

6.  Concluding  remarks  and  future  directions 

This  paper  reviewed  recent  development  on  MNF  for  heat 
transfer  applications.  Recent  studies  have  dominated  the  ther¬ 
momagnetic  convection  and  the  development  of  devices  by 
means  of  thermomagnetic  effects.  Some  researchers  initiated 
the  investigations  on  other  aspects  of  heat  transfer  such  as  the 
enhancement  and  control  of  thermal  conductivity  of  MNF  in  the 
presence  of  an  external  magnetic  field.  The  use  of  MNF  in  heat 
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transfer  applications  appears  promising.  However,  the  develop¬ 
ment  of  this  area  is  still  challenged  by  many  aspects  ranging 
from  the  preparation  and  characterization  of  MNF  to  practical 
applications  through  the  understanding  of  the  mechanisms 
responsible  for  the  observed  heat  transfer  enhancement.  In  order 
to  design  systems  with  enhanced  heat  transfer  and/or  heat 
transfer  effectiveness  (heat  transfer  coefficient)  using  a  MNF  as 
a  heat  transfer  fluid,  the  understanding  of  the  mechanisms  and 
the  physical  phenomena  is  needed.  The  problems  of  hydrody¬ 
namics  and  heat/mass  transfer  have  been  dealt  with.  However, 
some  problems  related  to  the  stability  of  MNF  in  the  presence  of 
magnetic  and  gravitational  fields,  heat  transfer  mechanisms  and 
the  influence  of  formulation  methods  on  the  thermophysical 
properties  are  still  needed  to  be  solved. 

Recent  studies  on  MNF  preparation  mainly  focused  on  the 
preparation  of  temperature  sensitive  super-paramagnetic  nano¬ 
particles  and  their  dispersion  in  appropriate  carrier  fluids.  In  an 
ideal  MNF,  each  particle  appears  superparamagnetic  with  the 
magnetic  moment  m.  However,  a  real  MNF,  as  used  in  some 
experiments,  is  a  poly-dispersed  system  and  may  contain  a 
small  fraction  of  single-  and/or  multi-domain  ferro-/ferri-mag- 
netic  nanoparticles.  The  presence  of  these  coarser  particles 
could  lead  to  a  poor  stability  of  MNF  and  have  severe  effects 
on  the  thermophysical  properties  of  the  fluid  and  the  heat 
transfer  effectiveness  of  systems  as  well.  In  addition,  the 
coarser  particles  in  the  MNF  samples  in  the  presence  of 
magnetic  field  lead  to  the  formation  of  long  and  irreversible 
chain-like  clusters,  which  have  detrimental  effects  on  the 
thermophysical  properties  (such  as  thermal  conductivity  and 
viscosity)  of  the  fluid.  A  higher  thermal  conductivity  could  be 
expected  for  MNF  prepared  with  the  smaller  particles.  It  was 
interesting  to  note  that  the  polydispersity  of  magnetic  nanopar¬ 
ticles  could  also  affect  directly  the  current  of  thermomagnetic 
convection.  Future  investigations  on  the  preparation  of  MNF  should 
focus  on  the  control  of  size/size  distribution  of  magnetic  nanoparti- 
cles.  Another  aspect  for  further  development  is  the  particle  morphol¬ 
ogy  control  during  the  preparation.  It  is  well  known  that  the 
magnetic  properties  of  magnetic  nanoparticles  and  hence  the  cluster 
formation  in  magnetic  nanofluids  could  be  affected  by  the  particle 
morphology.  The  control  of  particle  size/size  distribution  together 
with  the  particle  morphology  may  lead  to  the  preparation  of  well- 
engineered  MNF  with  desired  properties  and  long-term  stability  and 
further  on  the  development  of  devices  using  MNF  for  practical 
applications. 

Some  studies  have  shown  that  the  thermal  conductivity 
increased  with  the  particle  volume  fraction  in  the  absence  of 
magnetic  field.  It  is  clear  that  the  intensive  investigations  are 
needed  to  understand  the  mechanisms  and  the  influence  of 
other  parameters  (such  as  surfactant  layer,  particle  morphology, 
temperature,  etc.)  on  the  thermal  conductivity  of  MNF  in  zero 
magnetic  fields. 

The  investigations  on  thermal  conduction  in  the  presence  of 
magnetic  field  have  shown  a  promising  and  potential  application 
in  various  fields  such  as  smart  cooling.  The  anomalous  thermal 
conductivity  observed  was  explained  by  the  formation  of  chain¬ 
like  structures  in  MNF  under  the  magnetic  field.  The  under¬ 
standing  of  these  structures  and  their  manipulation  by  an 
external  magnetic  field  becomes  the  main  challenge  of 
the  future  investigation  on  thermal  conductivity  enhancement 
in  MNF.  The  control  of  the  self-assembly  processes  requires  the 
understanding  of  the  relationship  between  the  magnetic  proper¬ 
ties  and  other  particle  characteristics  such  as  particle  size/size 
distribution,  particle  morphology  as  well  as  the  applied  mag¬ 
netic  field.  The  future  investigations  on  thermal  conductivity 
enhancement  in  the  presence  of  magnetic  field  should  focus  on 
the  following  aspects: 


a)  understanding  of  the  interplay  among  the  parameters  such  as 
particle  size/size  distribution,  particle  morphology  and  coating 
shell  and  external  magnetic  field  as  well  as  its  influence  on 
particle  clustering  in  MNF; 

b)  development  of  thermal  conductivity  models  based  on  particle 
clustering,  clusters  morphology  and  distribution; 

c)  understanding  and  optimization  of  the  synergic  effect  of 
thermal  conductivity  and  viscosity  of  MNF  in  the  presence  of 
magnetic  field; 

d)  development  of  the  thermal  conductivity  and  viscosity  of  MNF 
with  magnetic  nanoparticles-carbon  nanotube  composites,  or 
other  novel  materials;  and 

e)  development  of  cooling  based  devices  and  the  potential 
application  of  the  thermal  conductivity  tuning  in  energy 
conversion  systems  and  heat  exchangers. 

The  thermomagnetic  convection  received  a  particular  atten¬ 
tion,  and  the  results  on  the  corresponding  numerical  simulations 
and  modeling  are  shown  in  the  literature.  Future  investigation 
should  focus  on  the  experimental  validation  of  the  existing 
models.  In  addition,  the  temperature-sensitive  magnetic  nano¬ 
fluid  where  the  thermomagnetic  convection  is  mainly  enhanced 
through  the  pyromagnetic  coefficient  is  investigated.  However, 
recent  studies  have  shown  that  the  thermomagnetic  convection 
could  be  strongly  affected  by  the  magnetoviscous  effect  as  well  as 
the  magnetic  field  effect.  In  future,  this  aspect  should  be  also 
investigated  for  both  temperature  sensitive  magnetic  nanofluids 
and  magnetite-based  nanofluids.  Another  aspect  for  further 
development  is  the  thermomagnetic  convection  affected  by 
modulated  magnetic  fields. 

It  was  demonstrated  that  the  solar  radiation  could  be  com¬ 
pletely  absorbed  in  a  magnetic  nanofluid  layer  with  about  10-mm 
thickness  [117].  This  offers  a  possibility  of  using  MNF  as  a  heat 
transfer  medium  in  solar  installations.  However,  studies  regarding 
the  parameter  effects  on  the  overall  heat  transfer  coefficient  and 
the  fluid  flow  within  the  system  still  are  needed.  Future  work 
should  focus  on  the  investigation  on  the  optical  and  thermal 
properties  of  MNF  as  well  as  their  possible  influence  on  the 
performance  of  solar  collectors,  especially  in  the  case  of  direct 
absorbing  solar  collectors.  Another  interesting  aspect  may  be  a 
possible  design  of  solar  thermosyphons  where  the  fluid  flow  and 
heat  transfer  process  may  be  completely  controlled  by  an  external 
magnetic  field. 

The  use  of  MNF  for  heat  transfer  enhancement  is  promising. 
This  field  has  made  a  significant  progress  for  some  challenging 
applications  such  as  smart  cooling,  automatic  energy  devices, 
automatic  cooling  devices,  etc.  However,  it  is  still  on  the  way  to 
go  before  some  novel  devices  using  MNF  as  a  heat  transfer 
medium  may  satisfy  the  practical  applications. 
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